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1. INTRODUCTION 

Theoretical calculations predict that cloud reflectance in near-infrared windows such as 
those at 1.6^/m and 2.2*/m should give lower reflectances than at visible wavelengths (Pollack et 
al„ 1978; Hansen and Pollack, 1970; Twomey, 1971). The reason for this difference is that ice and 
liquid water show significant absorption at these wavelengths, in contrast to the nearly 
conservative scattering at wavelengths shorter than l//m. In addition, because the amount of 
absorption scales with the path length of radiation through the particle, increasing cloud 
particle size should lead to decreasing reflectances at 1.6//m and 2.2//m. Measurements at these 
wavelengths to date, however, have often given unpredicted results. Twomey and Cocks (1982) found 
unexpectedly high absorption (factors of 3 to 5) in optically thick liquid water clouds. Curran 
and Wu (1982) found unexpectedly low absorption in optically thick high clouds, and postulated the 
existence of supercooled small water droplets in place of the expected large ice particles. We 
will examine the implications of the FIRE data for optically thin cirrus. 

2 RESULTS 

The Landsat satellite has spectral bands at 0.83*/m, 1.65//m, and 2.21//m which cover this 
range of variation in cloud absorption. Each pixel has a nominal spatial resolution of 28.5 
meters. Figure 1 gives the region covered by the Landsat data over Lake Michigan on October 28, 
1986. Figure 2 gives the Landsat measured nadir reflectance ratio R(2.21//m)/R(0.83^m) for the 
58.4 km square analysis region (solid line) in Fig. 1. At 15:38:30 UTC the King Air aircraft took 
a direct sample of the cloud particles on an oil covered slide. The sample is shown in Fig. 3b 
and is dominated by water droplets with a mean radius of about 4//m. This sample corresponds to a 
reflectance ratio of about 0.75 found in the Landsat data at location 1 in Fig. 2. There is a 
time difference of 15 minutes between the King Air and Landsat observations. The liquid water 
regions of this cloud, however, appear to have been colloidally stable (Heymsfield et al, 1989). 
A second direct cloud particle sample was collected at 15:52 UTC. This sample is shown in Fig. 3a 
and contains only ice particles (broken spatial plates and some columns, 20 to 300Am in length). 
This second sample corresponds to a reflectance ratio of about 0.4 found in the Landsat data at 
location "2" in Fig. 2. Note that the reflectances used to derive the image in Fig. 2 are not 
corrected for surface reflectance. In this case the reflectance ratios are a mixture of clear and 
cloudy signatures. 


rsU L* i * - 1 ^ t 


FLA?:X MOT FILMED 


369 



Given this qualitative agreement between the satellite and aircraft data, the next step is 
to test the quantitative agreement along the King Air aircraft track. The Landsat radiance data 
are spatially averaged to 1 km resolution, sampled every 0.5 km along the King Air groundtrack. 
Cloud reflectances are then corrected for surface reflectance effects as in Platt et al (1980). 
The 1.65//m and 2.21//m channels are found to require less than 0.01 correction for surface 
reflectance. The correction for the 0.83//m channel is less than 0.05. 

Figure 4 compares theoretical radiance calculations using the Finite Difference method 
(Suttles, 1981, 1985; Barkstrom, 1976) with the measured nadir cloud reflectance at 0.83//m and at 
both 1.65//m and 2.21//m along the King Air groundtrack. Calculations use a solar zenith angle of 
60 . The phase function for ice particles is taken from the laboratory measurements of 
Volkovitskiy et al (1980). The phase function for water particles is taken from theoretical Mie 
calculations with an effective radius of 3.8//m. Figures 4a and b give results for 1.65^m. 
Figures 4c and 4d give results for 2.21//m. The Landsat data are shown with symbols indicating the 
corresponding portion of the King Air track shown in Fig. 2. 

It is evident that there are two distinct populations of cloud particles along the 88 km 
track. The high reflectance ratio values in Fig. 2 (15:38:10-15:39:09 UTC and 15:49:40-15:50:39 
UTC) appear along the diagonal of nearly equal reflectance at the two wavelengths and are 
consistent with water droplets or ice spheres with radius less than 7.5//m. The remaining data 
indicate larger particles of about 60//m radius. Examination of the 1.65/ym versus 2.2\^m data 
given in Fig. 4 indicates that the large particles are ice. An assumption of liquid water for the 
large particles would give inconsistent particle sizes at 1.65 and 2.21//m. An assumption of ice 
gives consistent particle size in the two wavelengths. The small particles are too small to 
reliably distinguish ice from liquid water phase for these optically thin clouds. 

Figure 5 gives the King Air particle size distributions using the combined FSSP, 2D-C, and 
2D-P probes. F^r the 2D-C and 2D-P probes, particle size is calculated as a sphere with cross- 
section area n r equal to the area of the particle image in the 2-D probe. For compact non- 
spherical particles, this specification is similar to using equivalent volume spheres (Pollack and 
Cuzzi, 1980). Pollack and Cuzzi (1980) found that for large size parameter x = 2nr/A » 1, and 
moderate absorption 2n’x < 1) (where n’ is the imaginary index of refraction), equivalent volume 
spheres are most accurate for absorption efficiency determination. For the 1.65^/m and 2.2 l^m 
spectral bands, the appropriate radius range would be from about 3^m to 400//m. Given the compact 
particle habits observed in the microphysical data, the use of equivalent cross-section area 
spheres should be reasonably accurate. Large aspect ratio particles would lead to an overestimate 
of the true particle volume, and therefore an overestimate in particle absorption at 1.65 and 
2.21/sm. Such particles, however, were rarely noted in the data. 

Four characteristic size distributions are given in Fig. 5. The data for 15:38:10 to 

15:39:09 UTC cover the high ratio of R(2.21)/R(0.83) found in Fig. 2 near the location "1" in the 
figure. The microphysical data are dominated by small water droplets with radius about 4//m. The 
next section of the flight track (15:39:10 to 15:45:24 UTC) shows a peak at about 150//m, but no 
water droplets. The third section (15:49:40 to 15:50:39 UTC) has the smallest particle 
concentrations, and is taken from the high ratio of R(2.21)/R(0.83) found just before the end of 
the King Air track in the Landsat image. The Landsat data imply small particle sizes, while the 
aircraft finds no small drops in the FSSP probe. Spinhirne and Hart (1989) noted from the ER-2 
lidar data (ground track shown in Fig. 2) that the mixed phase cloud occurred in vertically thin 
layers (100 - 200 meters thick) at heights between 7.3 and 8.0 km. The lidar depolarization data 
at location in Fig. 2 verifies the existence of a mixed ice/liquid water phase cloud layer at 
7.3 km altitude (Spinhirne and Hart, 1989), the position of the King Air at 15:38:30 UTC. The 
King Air altitude at 15:50 UTC is 7.0 km, which is below the lidar detected altitudes for mixed 
phase cloud. It is likely that the King Air data at 15:50 missed the liquid water layer. We 
conclude that the aircraft microphysics and Landsat reflectances are in qualitative agreement, 
subject to uncertainties in the vertical variation of cloud microphysics and temporal evolution of 
the cloud field. 

The quantitative comparison of aircraft and radiometrically derived particle size requires 
the determination of an effective mean particle radius. Figure 6 gives the Landsat 2.21/0.83//m 
cloud reflectance ratio versus effective radius r . We define r = / r N(r)dr / f rN(r)dr, where 
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N(r) is the size distribution derived using the King Air microphysics data. This effective radius 
is a cross-section area weighted radius and has been shown useful in characterizing overall 
radiative properties of a particle size distribution. The number densities are averages over 
5-second intervals, which results in size distributions representative of 500-meter sections of 
the cloud field. 

Figure 6 indicates a significant but poor correlation between aircraft particle size and 
the 2.21/0. 83//m reflectance ratio. There appear to be three clusters of data, one with r g of 
about A pm, and two with r about 200pm. The apparently anomalous data with large particle size 
and large reflectance ratit? (15:50 UTC) are the liquid water layer missed by the King Air as 
discussed above. Recall that the data in Fig. 4 also clustered in two particle size groups, water 
droplets with radius less than 7.5 micron, and ice particles with radius about 60//m. While the 
water droplets appear consistent between the two data sources, the ice particles are in 
substantial disagreement, the radiative measurements indicating a smaller particle size by about a 
factor of 3. In order to understand this discrepancy, the uncertainties inherent in such a 
comparison are examined below. 

Uncertain Index of Refraction: factor of 2. 

Warren (1984) estimates that the uncertainty of the imaginary index of refraction for ice 
in the 1.4 to 2.%pm spectral region is a factor of 2. 
it. Uncertainty in the scattering phase function: - 30%. 

Use of the Mie scattering phase function in place of the VPP gave particle sizes about 60% 
smaller than those predicted in Fig. 4. This is an extreme test of the phase function. 

Hi. Uncertainty in the use of a single particle radius to represent an entire size 

distribution: * 15% at r < 200 pm. 

Mie calculations were fun to determine single scatter albedo as a function of r for the 
King Air size distributions (5 second averages) for each wavelength. These complete calculations 
were then compared to Mie calculations using a single particle size (as in Fig. 4). 

iv. Uncertainty in converting 2D image area to equivalent sphere: Unknown. 

The good agreement between the King Air median mass-weighted diameter determined as a 
function of crystal habit and 2 r determined using particle cross-section area (Fig. 6) argues 
that this error is of secondary importance. 

v. Inadequate microphysical sampling: Unknown. 

The size comparisons in Fig. 6 are given using the King Air measurements in the lower cloud 
layer at 7.3 km. Measurements of the upper cloud layer by the Sabreliner show mean particle 
radius decreasing from 200pm at 9 km to 40|/m at 11 km. Note that a 20pm radius is the smallest 
particle size measured by the 2D-C probe at the Sabreliner airspeed. Many small particles in the 
upper cloud layer are undoubtedly missed by the 2D-C and 2D-P probes (Heymsfield et al, 1989). In 
support of this concern, Spinhirne and Hart (1989) estimate a mode radius of 20pm for the upper 
cirrus layer on October 28 using integrated lidar backscatter and ll|/m emittance. For the lower 
cloud layer at 7-8 km. Fig. 5 shows that all of the 2-D size distributions at 7.3 km have maximum 
number density for the smallest particle size measured by the 2D-C probes. While the effect of 
these unsampled small ice particles on the determination of r g is unknown, it seems reasonable 
from the size distributions shown in Fig. 5, that this effect alone might explain the discrepancy 
between 60 and 200pm values for r . 

3. CONCLUSIONS 

The FIRE October 28, 1986, data provide a unique opportunity to compare measured and 
theoretical cloud properties for cirrus clouds. Overall impressions are: 

1. The lower cloud layer (7-8 km) appeared to dominate the cloud radiative properties as 
viewed by the ER-2 and Landsat. This result is consistent with the King Air and Sabreliner 
microphysical measurements and with ER-2 lidar observations. 

2. Particle size inferred using Landsat cloud reflectance at 0.83//m, 1.65//m and 2.2\pm gave 
good agreement with the King Air cloud particle samples for portions of the cloud field dominated 
by small water droplets with r = 3.8//m. For the larger ice crystals, the radiation measurements 
determined an r of about 60//m, compared to about 200//m determined using the King Air FSSP , 2D-C, 
and 2D-P proves. We conclude that the discrepancy is caused by two uncertainties. First, ice 
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particle sizes below about 20//m are not detected by the aircraft probes. The particle number 
densities are maximum at the smallest particles sensed by the 2-D probes, indicating the presence 
of smaller ice particles, even in the radiatively dominant lower cloud layer at 7-8 km. Second, 
uncertainties in the imaginary index of ice for the 1.65 and 2.2l#/m spectral bands causes an 
uncertainty of a factor of two in the Landsat derived particle size. The effect of small ice 
particles in the upper cirrus layer at 9-11.6 km is estimated to be of secondary importance. 

3. Recommendations for future cirrus experiments include improved measurement of ice particle 
concentrations for sizes between 1 and 50/;m, improved sampling of the vertical variation of cloud 
microphysics, more accurate radiometric calibration of aircraft radiometers, and more accurate 
values for the imaginary index of ice between 1.5 and 2.5#/m wavelengths. 
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Landsat 5: October 28, 1986 
15:53 UTC 
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LANDSAT, R(2.1 nm)/R(0.83 ym), 15:53:12 GMT 
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Fig. I. Location of the Landsat image area for the 
study. Solid line box gives the 58,4 km square area over 
Lake Michigan analyzed and shown in Fig. 2. 


Fig. 2. Landsat reflectance ratio, R(2.2l//m)/R(0.83>/m) 
over the analysis area. Aircraft tracks and observation 
times are also given in the figure. Data in Figs. 5 and 
6 is taken along the King Air ground track. 
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Landsat Particle Size Estimation 
along the King Air track 


King Air, 15:51:57 UTC. 7.0 km 
-27.8 C 



King Air, 15:38:30 UTC, 7.3 km 
-29.8 C 
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Fig, 3. Photographs of cirrus particles 
collected by the King Air on oil coated slides. 
3a shows ice crystals collected at 15:51:57 UTC 
and corresponds to location "2“ in Fig. 2. 3b 
shows water droplets collected at 15:38:30 UTC, 
and corresponds to location *1“ in Fig. 2. 


Fig. A. Measured and calculated cloud nadir reflectance. Landsat 
observations are taken along the King Air ground track seen in Fig. 2. 
Theoretical calculations use the VPP ice scattering phase function and 
Mie single scatter albedos as a function of particle radius. 4a and 4b 
give results for R(0.83j/m) vs. R().65jum) for liquid water (a) and ice 
(b) refractive index. 4c and 4d give results for R(0.83//m) vs. 
R(2.2lim) for liquid water (c) and ice (d) refractive index. 


Particle Size Distribution 



Fig. 5. King Air measured cirrus size distributions 
using the FSSP, 2D-C and 2D-P probes. Particle size for 
2D probes is that of a sphere with equivalent cross- 
section area to the particle 2-D image. Size 
distributions are averaged over the time intervals shown 
in the figure. 


Landsat Reflectance Ratio vs. 



Fig. 6. Comparison of Landsat measured cloud reflectance 
ratio R(2.2l*m)/R(0.83Am) with the King Air determined 
effective radius r^. 
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